June 1995 

IFUNAM FT 95-78, KL-TH 95/18 



Exact Solvability of the Calogero and 
Sutherland Models g 

by 

Werner Riihl" 

Fachbereich Physik - Theoretische Physik, University of Kaiserslautern, 
Kaiserslautern D-67653, Germany 

Alexander Turbiner 6 

Institute) de Fisica, UN AM, Apartado Postal 20-364, 01000 Mexico D.F., 

Mexico 

ABSTRACT 

Translationally invariant symmetric polynomials as coordinates 
for iV-body problems with identical particles are proposed. It is 
shown that in those coordinates the Calogero and Sutherland N- 
body Hamiltonians, after appropriate gauge transformations, can 
be presented as a quadratic polynomial in the generators of the 
algebra sIn in finite-dimensional degenerate representation. The 
exact solvability of these models follows from the existence of the 
infinite flag of such representation spaces, preserved by the above 
Hamiltonians. A connection with Jack polynomials is discussed. 
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The Calogero and Sutherland models are remarkable iV-body classical 
and quantum problems with a pair-wise interaction that are completely in- 
tegrable and at the same time are exactly solvable in their quantum versions 
[|T], (see also an excellent review ||). However, unlike the notion of inte- 
grability, which had been well established, for a long time the meaning of 
exact solvability remained intuitive, intrinsically folkloric. Just recently, an 
attempt to establish a certain definition of exact solvability was made [[|. 

Let us take an infinite set of finite-dimensional spaces Vo, Vi, V%, . . . with 
an explicit basis. Assume that they are embedded one into another, forming 
the infinite flag of spaces: Vo C V\ C V2 C — The simplest example of such 
a flag of spaces V n is given by a linear space of polynomials in one variable 
of degree not higher than n : 

V n = span{l, x, x 2 , . . . x n } . (1) 

Then one can give a definition saying that a linear operator h is named 
exactly solvable if it preserves the infinite flag of the spaces V n : 

h: V n ^V n , 71 = 0,1,2,... (2) 

Such a definition immediately leads to the conclusion that, in the basis where 
all spaces V n are naturally defined, the operator h in matrix form is given by 
a triangular matrix. This does not contradict h being self-adjoint. 

Now one can state that our task is to find the flags of different spaces V n , 
where the Hamilton operators act. Once we succeed in this task, we end up 
with the exactly solvable Hamilton operators. 

There are different ways to choose the spaces V n . One of the most natural 
choices is to identify the space V n with a finite-dimensional representation 
space of a certain semi-simple Lie algebra. Then it is evident that an exactly 
solvable operator h generically should be a nonlinear combination of gener- 
ators of the Cartan and negative root generators. For instance, the space 
(1) corresponds to the (n + l)-dimensional representation of the algebra SI2 
realized by first-order differential operators acting on the real line. 

Of course, such a choice of V n is very particular. However, it is quite 
amazing that all known one- dimensional exactly solvable quantum mechani- 
cal problems such as the harmonic oscillator, the Coulomb problem, etc., can 
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be interpreted in such a manner as well as all exactly solvable discrete oper- 
ators having Hahn, Kravchuk, Pollachek, Charlier and Meixner polynomials 
as their eigenfunctions ||. The purpose of the present paper is to show that 
the iV-body Calogero and Sutherland models also belong to this type being 
matched to finite-dimensional representions of the algebra sl^. Thus the list 
of finite-dimensional exactly solvable problems is exhausted and we arrived 
at the conclusion that all of them are of the same type, being presented by 
a quadratic combination of the generators of a certain Lie algebra. 

1. The Calogero model 

The Calogero model fl[] is a quantum mechanical system of N particles on 
a line interacting via a pair-wise potential; it is defined by the Hamiltonian 

i N N 

Hc al = lEH + - 2 4+T,7^y 2 , (3) 

^ i=l j<i V x « x j/ 

where di = oj is the harmonic oscillator frequency and hereafter normal- 
ized to uj = 1, and g = v(y — 1) is the coupling constant. 

As was found by Calogero this model is completely integrable and 
possesses some very remarkable properties. Firstly, in order to be normaliz- 
able, all eigenfunctions must (up to a certain function as a common factor) 
be either totally symmetric or totally antisymmetric. Secondly, it turned 
out that the energy spectrum is that of N bosons or fermions interacting via 
harmonic forces only, but with a total energy shift proportional to v (for a 
review see ||). Furthermore, the multiplicity of degeneracies of the states is 
the same as in the pure oscillator problem without repulsion. It is also worth 
noting that for N = 3 (and for this case only) the problem can be solved by 
a separation of variables. 

Thus the problem (|3|) is not only completely integrable, but also exactly 
solvable: the spectrum can be found explicitly, in a closed analytic form. 

All eigenfunctions have a form (3 V e ~P(x), where (3 is the Vandermonde 
determinant, X 2 = Y^iLi x % \ v * equals either v or (1 — u), and P(x) is a com- 
pletely symmetric polynomial under the permutation of any two coordinates. 
Thirdly, after extracting the center-of-mass motion, the remaining operator 
is translationally invariant under: Xi — > Xi + a. 
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Now let us make a rotation of the Hamiltonian (|3|) || : 

h = -2(3~ v * HcaiF* = 

N N N 1 

i=l i=l j^i Xi X 3 

- N - u*N(N - 1) , (4) 

where g = u*(v* — 1), thereafter the constant term in (4) will be omitted. 
Introduce the center-of-mass 

N 



y = E 



X 



i=i 

and translationally invariant Jacobi coordinates J7[ (see also 0]) : 

1 ^ 

?/i = a;. - -77 E x i ' i = l,2,...,iV (5) 

fulfilling the constraint Y^=i Hi = 0. 

It is known that the appropriate eigenfunctions of the operator (4) are 
completely symmetric under the permutations of x-coordinates and, con- 
sequently, of ^-coordinates. Now, instead of the explicit symmetrization of 
the eigenfunctions of (4), we encode this feature in a coordinate system con- 
sisting of functions in the {xi} and {yi}, which are completely symmetric un- 
der the permutations. The eigenfunctions, depending on these coordinates, 
are then automatically symmetric. 

There exist many different types of symmetric polynomials. For our pur- 
pose, the most adequate set of symmetric polynomials is given by the ele- 
mentary symmetric (ES) polynomials ||: 

a i( x ) = E X i ' 

Ki<N 



°2(a0 = E 



l<i<j<N 

a s( x ) = E x i x j x k 
l<i<j<k<N 



<T N (x) = X!X 2 X 3 ...X N , (6) 
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where x G R. We notice that <J\{x) coincides with the center-of-mass coor- 
dinate Y. Those polynomials can be defined also by a generating function 



N N 

n(i+^) = 5>n(z)* n 

i=l n=0 



(7) 



with (To = 1. The ES functions can be used as coordinates on the subset 
EcR N 

< x 1 < x 2 < x 3 . . . < Xn < oo 

of R^. In turn, differential operators acting on E can be rewritten in terms 
of ES coordinates. If these differential operators are themselves symmetric, 
they can be extended from E to R^. 

It is worth emphasizing that the coordinates (6) allow us to avoid the 
problem of over-completeness of the basis, which appears if the Newton poly- 
nomials 



N 



s k = J2 



X: 



are used as the coordinates [|7| . 

Moreover, the Jacobian in the volume form 

dai A da 2 A ... A daw = JnAxi A dx 2 A . . . A dx 



N 



(9) 



is 



N 



x l7 x 2 , ... ,x N ) (10) 

with f3 the Vandermonde determinant as before and [a] means the integer 
part of a. In turn we can express J N in terms of the polynomials Sk © 



J 



N 

S 3 



s l S 2 S3 

s 2 s 3 

S3 



SN-1 

sn 



sn-i sn 



S2N-2 
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These Sk (including those with k > N) can be expressed by the a n with the 
help of the generating function (0) 

00 ( — ] \n+l ( N \ 

Y, [ -^Sn(x)t n = \0g(J2vn(x)t n ) . (12) 
n=l H \n=0 / 

The coordinates (6) are not translationally invariant, but this property 
can also be included by considering the ES polynomials of the arguments (5), 

r n (x) = a n (y(x)) , n = 2,3, . . . , N , (13) 

as the relative coordinates and with the remaining <Ji(x) = Y as the coor- 
dinate of the center of mass. Certainly, once n > N, all r n (x) = 0. So r n 
implicitly contains the information on the value of N. 

Making quite sophisticated but straightforward calculations, one can get 
an explicit expression for the Laplacian in r-coordinates (7): 

N Q2 Q2 N Q2 N q 

4s S»? = w s? + ^ + S\' (14) 

where 

(N-j + l)(k-l) ^ . n/A 

Ajk = Tj-iTk-i + 2^ \k — J — 2£)r j+£ _ir fc _ £ _i 

£>max(l,fc~ i j) 

D {N -i + 2)(N -2 + 1) 

* = N Tl - 2 • 

Here we put r = 1, T\ = and r p = 0, if p < and p > N. 

Finally, after the extraction of the center-of-mass motion and omitting 
the constant terms, the operator (4) reduces to h re i for the relative motion 

N f)2 N 53 / 1 \ iV Q 

hr*= E ^^-2E^^-(^ + ^)E(^+2)(iV-z+l)r i _ 2 - 

(15) 



Now one can pose a question: would it be possible to rewrite (J15|) as 
an element of the universal enveloping algebra of a certain Lie algebra in a 
finite-dimensional representation? 
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Let us take the algebra gl^IL). One of the simplest representations of 
this algebra in terms of first-order differential operators is the following : 

Jr = J-, i = 2,3,...,JV, 
_ Q 

J i j = tiJj = t i 7^— , z,j = 2,3,...,iV, 

N Q 

p=2 UL P 

J+ = UJ\ i = 2,3,...,JV. (16) 

which acts on functions of t G R N_1 . One of the generators, namely J° + 
Y,p=2 Jp, P i s proportional to a constant and, if it is extracted, we end up 
with the algebra sZjv(R). The generators Jj^ form the algebra of the vector 
fields sZjv-i(R). The parameter n in ( |I6"D can be any real number. If n is a 
non-negative integer, the representation ([16]) becomes the finite-dimensional 
representation acting on the space of polynomials 

V n (t) = span{t" 2 ^C . . . t n N N : < £ m < n} . (17) 

This representation corresponds to a Young tableau of one row and n blocks 
and is irreducible. 

It is easy to see that the operator h re i for the relative motion can be 
rewritten in terms of the generators (|T6|). The representation of h re i in terms 
of generators flTED is 

h - V ) (N-j + l)(j-l) ( 2 _ ^ o o 

«re/ - 2^ \ N \ J j-i,j) z 2^ lJ j+e-i,j J j-e-i,j 




=i 

\ , (N-j + l)(k-l) o 

_r Z 2^ 1 J j-l,j J k-l,k 

N 



fc-1 ~> V 

/Xj — k + 2£) Jj +e _ 1 jJ k _ e _ l k > — 2 fcJ fc fc 
/ 1 \ * 

+ u *j 12( N ~k + 2){N-k + l)J° k _ 2>k (18) 
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where we identify Jq^ with J k and put J° fc equal to zero. Moreover the U 
in ( JL6D are replaced by the Tj. 

Such a rewriting ([Tj|) can be performed when the parameter n takes any 
value. Hence, the operator h re i possesses infinitely many finite-dimensional 
invariant subspaces V n (t), n = 0, 1, 2, . . . , and, correspondingly, preserves an 
infinite flag of the spaces V n (t). Therefore this operator is exactly-solvable 
according to the definition given above. The case N = 3 confirms the general 
hypothesis stated in Ref. || about the non-existence of exactly and quasi- 
exactly solvable problems in two-dimensional flat space without separability 
of variables. Also the operator (18) at iV = 3 is not contained in the lists of 
the (quasi)-exactly-solvable operators in R 2 presented in the papers 0. 

2. The Sutherland model 

The Sutherland model |2j (for a review, see ||) is a quantum mechanical 
system of iV particles on a line interacting via a pair-wise potential and 
defined by the Hamiltonian (cf. (^)) {di = ^-) : 

1 N a 
H Suth = - 7 :T,dl+ E • 2i, _ \ (19) 

Z i=l l<i<j<N Sm 2V 1 ' X 3J 

which is defined on the Hilbert space of functions over the torus (S\) xN . As 
in (3), g = v[y — 1) here is the coupling constant. 

Similarly to the Calogero model, this model is completely integrable and 
exactly-solvable - the spectrum can be found explicitly, in a closed analytic 
form. All eigenfunctions, in order to be normalizable, must (up to a certain 
function as a common factor) be either totally symmetric or totally antisym- 
metric and then have a form f3 u P(x), where 



p= n 

l<i<j<N 



. 1, 

sin - {Xi — r 



is a natural modification of the Vandermonde determinant, u* equals either 
v or 1 — v, while P(x) is a completely symmetric polynomial in coordinates 
z = exp ix under the permutation of any two coordinates. These polynomials 
P(z) are named Jack polynomials ||. It is evident that the Hamiltonian is 
translationally invariant under Xi —>■ Xi + a. 
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Now we make the gauge transformation of the Hamiltonian (|i9|) (see ||) 
h = -2p- v *H Suth (3 v * 

N 

= d\ + iv* cot |(xj — Xj)(di — dj) + const (20) 

i=l l<i<j<N 

with the same relations between g, z/, z/ 1 as in the Calogero model. The op- 
erator h has Jack polynomials as eigenf unctions. 

In order to realize a possible hidden algebraic structure, one should in- 
troduce a torus analogue of the symmetric polynomial coordinates (7). Let 
us consider the torus 

Zir = (£i) xJV 

with standard coordinates 

Zl (x) = e ix \ i = l,...N, (21) 

< Xi < 2tt . 
Introduce the ES polynomials (see (6)) 

?7„(x) = a n (w) , (22) 

where 

^ = e ty '(x), Vi = Xi- jy<ri{x) ■ (23) 

Obviously the coordinates (|22|) are translationally invariant under — > £j + a 
and 

= - „ • (24) 

6v(z) 



Those variables possess a remarkable symmetry 

<0*0 = 77iV-n(z) (25) 

what implies that for even N the variable r]N/2(x) is real. 
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Finally, complete system of coordinates we take 

{£ N ,VnA<n<(N-l)} , (26) 

where £n describes the motion of the center-of-mass. In these coordinates 
the Laplacian has the form 

-A = nU n -£-) + £ B jk - 



with 



d£ N J drjjdrjk 

+ Vp^-^ (27) 

Bjk = ^ VjVk + ( k ~j- tyVj+m-i (28) 

Z>max(l,fc— j) 

(cf. ([14])) and here rjo — T]n — 1- The center-of-mass motion can be separated 
and after omitting the constant, additive term the Hamiltonian of the relative 
motion takes the form 

N-l Q2 
— Kel = J! Bjk 



dr]jdr] k 



j,k=l "'13 

+ (-*4)|V-0^. (29) 
With tj replaced by rjj-i in (|16"D we obtain 

- Kei = | ^ ~ ( J i.i ^ 2 ~ 2 51 / J i+' J J i-' J 

l<fc<j<iV-l I JV z=i J 

JV-l 

+ ^/(iV-Z)^ (30) 
z=i 

Here we identified Jp fc and J% k with J^. Like (18), the operator (30) at 
N = 3 is not contained in the lists of the (quasi)-exactly-solvable operators 
in R 2 presented in the papers M. 
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In conclusion, it is worth emphasizing that the Laplace operators (18), 
(27) (from the physical point of view it can be treated such as we study a free 
motion without an interparticle interaction), after extraction of the center- 
of-mass motion, are represented as a quadratic combination in the generators 
(|16D of the algebra sIn(R). The operator (29) realizes a new representation of 
Jack polynomials in the coordinates (22). This representation is much simpler 
than the standard representation (20) allowing easy an explicit construction 
of the first Jack polynomials. 
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